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FOREWORD 


Methodology  and  Characterization  Studies  during  fiscal  years 
1977  and  1978  included  gravity  surveys  in  ten  valleys  in  Arizona 
(five)  ,  Nevada  (two) ,  New  Mexico  (two) ,  and  California  (one) . 
The  gravity  data  were  obtained  for  the  purpose  of  estimating  the 
gross  structure  and  shape  of  the  basins  and  the  thickness  of  the 
valley  fill.  There  was  also  the  possibility  of  detecting 
shallow  rock  in  areas  between  boring  locations.  Generalized 
interpretations  from  these  surveys  were  included  in  Fugro 
National's  Characterization  Reports  (FN-TR-26a  through  e)  . 

During  the  FY  77  surveys,  the  measurements  were  made  to  form  an 
approximate  one-mile  grid  over  the  study  areas,  and  contour  maps 
showing  interpreted  depth  to  bedrock  were  made.  In  FY  79,  the 
decision  was  made  to  concentrate  the  available  funds  on  the 
basic  Verification  Program  to  verify  and  refine  suitable  area 
boundaries.  This  decision  resulted  in  a  reduction  in  the 
gravity  program.  Instead  of  obtaining  gravity  data  on  a  grid, 
the  reduced  program  consisted  of  obtaining  gravity  measure¬ 
ments  along  profiles  across  the  valleys  where  Verification 
Studies  were  also  performed. 

The  Defense  Mapping  Agency  (DMA),  St.  Louis,  was  also  requested 
to  provide  gravity  data  from  their  library  to  supplement  the 
gravity  profiles.  For  Big  Smoky,  Reveille,  and  Railroad  val¬ 
leys,  a  sufficient  density  of  library  data  is  available  to 
permit  construction  of  interpreted  contour  maps  instead  of 
two-dimensional  cross  sections. 

In  late  summer  of  FY  79,  supplementary  funds  became  available  to 
begin  data  reduction.  At  this  time,  inner  zone  terrain  correc¬ 
tions  began  on  the  library  data  and  the  profiles  from  Big  Smoky 
Valley,  Nevada,  and  Butler  and  La  Posa  valleys,  Arizona.  The 
profile  data  from  Whirlwind,  Hamlin,  Snake  East,  White  River  and 
Garden  Coal  valleys,  Nevada  were  available  from  the  field  in 
early  October,  1979. 

A  continuation  of  gravity  interpretations  has  been  incorporated 
into  the  FY  80  contract  and  the  results  are  being  summarized  in 
a  series  of  valley  reports.  The  reports  covering  Nevada-Utah 
gravity  studies  will  be  numbered,  "FN-TR-33-”,  followed  by  the 
abbreviation  for  the  subject  valley.  In  addition,  more  detailed 
reports  of  the  results  of  FY  77  surveys  in  Dry  Lake  and  Ralston 
valleys,  Nevada  are  being  prepared.  Verification  Studies  are 
continuing  in  FY  80  and  gravity  studies  are  included  in  the 
program.  DMA  will  continue  to  obtain  the  field  measurements  and 
it  is  planned  to  return  to  the  grid  pattern.  The  interpretation 
of  the  grid  data  will  allow  the  production  of  contour  maps  which 
will  be  valuable  in  the  deep  basin  structural  analysis  needed 
for  computer  modeling  in  the  Water  Resources  Program.  The 
gravity  interpretations  will  also  be  useful  in  the  Nuclear 
Hardness  and  Survivability  (NH&S)  evaluations. 
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The  basic  decisions  governing  the  gravity  program  are  made  by 
BMO  following  consultation  with  TRW  Inc.,  Fugro  National  and 
the  (DMA).  Conduct  of  the  gravity  studies  is  a  joint  effort 
between  DMA  and  Fugro  National.  The  field  work,  including 
planning,  logistics,  surveying,  and  meter  operation  is  done  by 
the  Defense  Mapping  Agency  Hydrographic/Topographic  Center 
(DMAHTC) ,  headquartered  in  Cheyenne,  Wyoming.  DMAHTC  reduces 
the  data  to  Simple  Bouguer  Anomaly  (see  Section  A1.4,  Appendix 
A1.0).  The  Defense  Mapping  Agency  Aerospace  Center  (DMAAC) ,  St. 
Louis,  calculates  outer  zone  terrain  corrections. 

Fugro  National  provides  DMA  with  schedules  showing  the  valleys 
with  the  highest  priorities.  Fugro  National  also  recommended 
locations  for  the  profiles  in  the  FY  79  studies  within  the 
constraints  that  they  should  follow  existing  roads  or  trails. 
Any  required  inner  zone  terrain  corrections  are  calculated  by 
Fugro  National  prior  to  making  geologic  interpretations. 
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1.0  INTRODUCTION 

1.1  OBJECTIVE 

Gravity  measurements  were  made  in  Whirlwind  Valley  for  the 
purpose  of  estimating  the  overall  shape  of  the  structural  basin 
and  the  thickness  of  alluvial  fill  in  the  basin.  The  esti¬ 
mates  will  be  useful  in  modeling  the  dynamic  response  of  the 
basin  to  explosive  or  earthquake-generated  ground  motion 
and  in  evaluating  ground-water  resources. 

1.2  LOCATION 

Whirlwind  Valley  is  located  in  western  Millard  and  western  Juab 
counties,  Utah  about  30  miles  (48  km)  west  of  the  town  of  Delta 
(Figure  1).  The  valley  is  bounded  on  the  east  by  the  Little 
Drum  Mountains  and  on  the  west  by  the  House  Range  (Figure  2)  . 
It  culminates  northward  in  Swasey  Bottom  near  the  Juab  County 
line,  and  is  open  southeastward  toward  Lake  Sevier  and  the 
Sevier  Desert.  U.S.  Highway  6/50  crosses  the  southern  end  of 
the  valley. 

1 . 3  SCOPE  OF  STUDY 

The  Defense  Mapping  Agency  Hydrographic/Topographic  Center 
(DMAHTC)  made  gravity  measurements  along  five  cross-valley 
profiles  as  illustrated  in  Figure  2.  The  distance  between  the 
profiles  ranged  from  5  to  8  miles  (8  to  13  km)  between  the 
profiles.  The  sampling  interval  was  1  mile  (1.6  km)  over  the 
central  valley  and  1/4  mile  (0.4  km)  near  the  valley  margins. 
The  more  dense  sampling  was  used  near  the  valley  margins  to 
define  any  gradients  associated  with  boundary  faults  and  to 
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identify  short-wavelegth  anomalies  which  might  occur  if  bedrock 
was  shallow  along  the  margins  of  the  valley.  The  station 
locations  are  shown  in  Figure  3  and  the  station  data  are 
listed  in  Appendix  A2.0.  The  tolerance  for  establishing  station 
elevations  was  5  feet  (1.5  m) .  This  tolerance  for  elevation 
control  limited  the  gravity  precision  to  0.3  milligals. 
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WW031Q  Gravity  station  number 

OOO  Gravity  station  locatis 

WW-1  Prof  i  I  e  name 
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Area  suitable  for  hybrid  trench  and  vertical  shelter 
basing  modes.  Depth  to  rock  and  water  greater  than 
150  feet  (46m). 

Area  suitable  for  hybrid  trench  and- not  suitable  for 
vertical  shelter.  Depth  to  rock  and  water  greater 
than  50  feet  ( ?5m)  and  less  than  150  feet  (46m). 

Area  unsuitable  for  both  hybrid  trench  and  vertical 
shelter  basing  modes  as  determined  from  application 
of  depth  to  rock  and  water,  topographic/terrain,  and 
cultural  exclusions.  (See  Section  A2.0  in  Appendix 
for  details  of  exclusion  criteria.) 


Indicates  areas  of  exposed  rock. 


Contact  between  Tock  and  basin  fill. 


Whirlwind  verification  site  boundary,  FY  79 
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2.0  GRAVITY  DATA  REDUCTION 

DMAHTC  obtained  the  basic  observations  for  each  station  and 
reduced  them  to  Simple  Bouguer  Anomalies  (SBA)  as  described  in 
Appendix  A1.0.  Up  to  three  levels  of  terrain  corrections  were 
applied  to  convert  the  SBA  to  the  Complete  Bouguer  Anomaly 
(CBA).  First,  the  Defense  Mapping  Agency  Aerospace  Center 
(DMAAC)  used  its  library  of  digitized  terrain  data  to  calculate 
corrections  for  terrain  to  104  miles  (167  km)  from  each  station. 
The  computation  method  used  has  limitations  for  correcting 
terrain  effects  near  the  stations.  This  made  it  necessary  to 
apply  a  second  level  of  correction  for  some  stations.  For 
these,  a  ring  template  was  used  to  calculate  the  effect  of 
terrain  within  approximately  3000  feet  (914  m)  of  the  station. 
A  third  level  of  terrain  correction  was  applied  to  those  sta¬ 
tions  where  10  feet  or  more  of  relief  were  observed  within 
130  feet  (40  m)  of  the  station.  For  these  stations,  the  ele¬ 
vation  differences  at  a  distance  of  130  feet  along  six  direc¬ 
tions  from  the  station  were  measured  in  the  field.  These  data 
were  used  to  calculate  the  effect  of  the  very  near  relief. 
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3.0  GEOLOGIC  SUMMARY 

Whirlwind  Valley  lies  within  the  Basin  and  Range  physiographic 
province.  Rocks  forming  the  Little  Drum  Mountains  on  the 
eastern  side  of  the  valley  consist  of  early  Tertiary  volcanic 
flows  (Hintze,  1963;  Stokes,  1963).  The  House  Range,  which 
forms  the  western  boundary,  contains  faulted  eastward-dipping 
limestones  and  dolomites  with  interbedded  shales  of  Cambrian 
age.  Ordovician  age  limestones  with  interbedded  conglomerates 
crop  out  south  of  the  valley  in  the  southern  House  Range  west  of 
Sevier  Lake.  Tertiary  age  volcanic  rocks  similar  to  those  in 
the  southern  Little  Drum  Mountains  are  found  along  the  valley 
axis  in  Red  Knolls  and  Long  Ridge.  Tertiary  age  conglomerates 
overlie  these  volcanic  rocks  and  are  the  principal  units  exposed 
in  Long  Ridge. 

Hintze  (1963)  and  Stokes  (1963)  suggest  only  one  fault  in 
Whirlwind  Valley.  They  show  it  concealed  by  alluvium,  trending 
southeastward  into  the  north-central  part  of  the  valley.  The 
House  Range  to  the  west  is  considered  by  Gehman  (1958)  to  be  the 
gently  eastward  dipping  limb  of  a  faulted  anticline.  If  this  is 
the  case,  the  carbonate  rocks  exposed  in  the  House  Range  extend 
beneath  the  alluvium  in  Whirlwind  Valley,  and  presumably  under¬ 
lie  the  volcanic  rocks  in  the  Little  Drum  Mountains.  The  steep 
western  face  of  the  House  Range,  bounding  the  next  valley, 
represents  the  eroded  fault  scarp  along  which  uplift  has  taken 
place.  Few  faults  are  mapped  within  the  Little  Drum  Mountains 
although  relations  are  complex  because  of  the  variety  of  vol¬ 
canic  rock  types  present. 
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Surficial  basin-fill  deposits  are  described  in  the  Verification 
Studies  (FY  79,  FN-TR-27-1A) .  They  are  predominantly  alluvial 
fan  and  lacustrine  deposits.  The  lacustrine  deposits  are 
associated  with  Pleistocene  age  Lake  Bonneville  and  occur 
principally  in  the  southern  and  central  parts  of  the  valley. 
The  alluvial  deposits  consist  primarily  of  silty  sand  with 
lesser  amounts  of  gravelly  sand,  silt,  and  clay.  Younger 
alluvial  fans  are  generally  uncemented;  intermediate  alluvial 
fans  are  weakly  cemented;  and  older  alluvial  fans  are  moderately 
to  strongly  cemented. 

Published  values  for  densities  of  soil  and  alluvium  generally 
range  between  1.3  and  2.3  g/cm3,  but  sometimes  attain  greater 
values  under  compaction  (Grant  and  West,  1965).  Limestone 
densities  generally  range  between  2.3  and  2.8  g/cm3  with  the 
older  formations  having  the  greater  densities. 
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4.0  INTERPRETATION 

A  valley  filled  with  relatively  low-density  alluvium  will  create 
a  negative  gravity  anomaly.  Thus,  a  graph  of  gravity  across  a 
valley  is  often  U-shaped,  low  in  the  middle  of  the  valley  where 
the  fill  is  thickest,  and  high  on  the  ends  where  the  fill  thins 
and  disappears.  The  basis  for  interpretation  is  the  CBA  pro¬ 
file.  The  CBA  for  the  five  profiles  across  Whirlwind  Valley  are 
shown  in  the  top  portion  of  Figures  4  through  8. 

4.1  REGIONAL  -  RESIDUAL  SEPARATION 

The  CBA  contains  gravitational  field  components  contributed  by 
geologic  conditions  which  are  unrelated  to  basin  fill.  These 
are  known  as  regional  effects  and  must  be  removed  from  the  CBA 
to  obtain  the  gravity  contributions  made  by  the  valley  fill. 
A  regional  field  was  established  by  linear  interpolation 
between  the  CBA  values  at  bedrock  stations  on  opposite  ends  of 
the  profiles.  Where  only  one  end  of  a  profile  was  on  bedrock, 
the  regional  value  was  assumed  to  be  constant  across  the  valley. 
This  method  does  not  result  in  true  residual  values  (that  is, 
some  regional  effects  remain)  but  the  error  is  probably  small 
compared  to  the  large  residual  anomaly  values  produced  by  the 
valley  fill. 

The  regional  trend  used  for  each  profile  is  shown  together  with 
the  CBA  in  the  top  portion  of  Figures  4  through  8.  The  differ¬ 
ence  between  the  regional  and  the  CBA  is  the  residual  gravity 
anomaly  which  is  attributed  to  the  presence  of  the  alluvial  fill 
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in  the  valley.  The  residual  for  each  profile  is  shown  by  the 
crosses  (x)  in  the  center  portion  of  Figures  4  through  8. 

4.2  DENSITY  SELECTION 

The  construction  of  a  geologic  model  to  account  for  the  residual 
anomaly  requires  selection  of  density  values  representative  of 
the  basin  fill  and  of  the  underlying  basement  rock.  Detailed 
density  values  are  not  available.  Consequently,  generalized 
values  were  chosen  and  used  for  modeling  all  five  of  the  pro¬ 
files. 

The  density  value  selected  to  represent  basin  floor  was 
2.8  g/cm^.  This  value  was  chosen  because  the  Paleozoic 
carbonates  in  the  Great  Basin  are  generally  reported  to  be 
relatively  dense.  It  is  assumed  that  the  carbonate  rocks 
forming  the  House  Range  extend  across  Whirlwind  Valley  beneath 
the  alluvium. 

Average  in  situ  densities  for  basin  fill  from  four  borings 
sampled  at  20-foot  (6  m)  intervals,  between  100  and  160  feet 
(30  to  49  m)  ,  ranged  from  2.0  to  2.3  g/cm^.  The  density  value 
selected  to  represent  the  basin  fill  was  2.3  g/cm^.  This 
value  is  near  the  high  end  of  the  densities  measured  for  samples 
from  160-foot  borings  drilled  during  the  Verification  Studies  of 
Whirlwind  Valley  (FY  79,  FN-TR-27-I I ) .  The  high  end,  instead  of 
the  average,  of  the  measured  values  was  used  because  all  the 
basin  fill  samples  came  from  relatively  shallow  depths  and 
densities  of  alluvial  materials  usually  increase  with  depth.  If 
compaction  at  depth  has  caused  the  overall  effective  density  to 
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be  greater  than  2.3  g/cm^.  The  thicknesses  of  alluvium  given 
by  the  models  in  Figures  4  through  8  are  less  than  the  actual 
thicknesses. 

4.3  MODELING 

Modeling  of  the  cross-sections  of  the  basin  beneath  the  gravity 
profiles  was  done  using  a  two-dimensional  computer-modeling 
program.  The  model  chosen  for  this  analysis  appears  in  cross- 
section  as  a  set  of  1-km-wide  blocks  whose  top  is  at  surface 
elevation  and  whose  bottom  represents  the  alluvium-bedrock 
boundary.  The  elevations  at  the  bottoms  of  the  blocks  were 
adjusted  by  iterative  computation  until  the  computed  gravity 
anomaly  for  the  model  differed  by  less  than  1/2  milligal 
from  the  observed  residual  gravity  anomaly.  The  computed 
gravity  anomaly  from  the  final  model  is  shown  as  a  continuous 
line  in  the  center  portion  of  Figures  4  through  8.  The  result¬ 
ing  depth-to-bed rock  models  are  shown  in  the  lowest  section  of 
Figures  4  through  8. 

The  shapes  of  the  basement  surfaces  derived  from  this  modeling 
have  a  basin-like  character  with  few  complications.  Profile  1, 
on  the  north  (Figure  4) ,  however  does  show  some  basement  irre¬ 
gularities  which  may  indicate  fault  blocks  or  intrusives.  The 
concealed  fault  suggested  by  Stokes  (1963)  ends  at  about  the 
location  of  Profile  1.  Although  presumably  the  downthrown  side 
would  be  on  the  east,  the  gravity  profile  shows  thinner  fill  on 
the  east  rather  than  thicker.  The  figures  have  a  vertical 
exaggeration  of  125/12  times  (10.4);  therefore,  the  gentle 
slopes  appear  steep. 


FN-TR-3  3-WW 


12 


Proceeding  southward,  Profile  2  (Figure  5)  shows  a  smooth  basin 
shape.  Profile  3  (Figure  6)  shows  half  a  basin.  Red  Knolls, 
which  lies  at  the  eastern  end  of  Profile  3,  has  little  or  no 
gravity  expression.  Consequently,  it  is  interpreted  to  be 
a  thin  layer  of  volcanic  flow  material  overlying  the  basin  fill. 
Profile  4  (Figure  7)  ,  also  shows  the  basin  floor  to  be  deepest 
under  the  western  Red  Knolls,  becoming  shallower  toward  the 
south  end  of  the  knolls.  The  gravity  profile  also  appears  to  be 
unaffected  by  Long  Ridge,  two  miles  to  the  south.  These  two 
profiles  indicate  that  neither  the  Red  Knolls  volcanics  nor 
the  Long  Ridge  conglomerates  are  significantly  more  dense  than 
the  other  basin  fill,  and  that  they  do  not  lay  directly  on 
bedrock.  Profile  5  (Figure  8),  which  passes  within  a  mile  of 
Sevier  Lake,  shows  a  bedrock  high  between  the  lake  and  Long 
Ridge.  The  calculated  bedrock  depth  is  about  240  feet  (73  m)  . 
This  high  point  may  be  a  structural  divide  between  Whirlwind 
Valley  and  the  Sevier  Desert  Basin. 

Figures  4  through  8  show  that  the  valley  fill  thickens  by  about 
400  feet  to  500  feet  (122  to  152  m)  per  kilometer  outward  from 
bedrock.  The  greatest  depth  of  fill,  about  2000  feet  (610  m) , 
is  beneath  Profile  1.  However  a  more  typical  depth  is  about 
1000  feet  (305  m) .  These  calculated  depths  are  very  dependent 
upon  the  choice  of  densities;  for  example,  a  one  percent 
increase  of  assumed  fill  density  for  the  model  makes  a  five 
percent  decrease  in  the  density  contrast  and,  therefore,  a 
five  percent  increase  in  the  calculated  thickness  of  fill. 
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5.0  CONCLUSIONS 

Whirlwind  Valley  is  interpreted  to  be  a  relatively  shallow 
structural  valley.  It  has  the  long,  narrow  shape  typically 
created  by  range  bounding  faults  in  the  Basin  and  Range  Physio¬ 
graphic  province. 

However,  if  such  faults  do  exist  in  Whirlwind  Valley,  they  must 
be  relatively  small  because,  they  are  not  visible  at  the  surface 
and  are  not  clearly  defined  by  the  gravity  profiles. 

The  interpretations  of  the  two  profiles  (3  and  4)  near  Red 
Knolls  and  Long  Ridge  suggest  that  these  topographic  features 
are  not  basin  boundaries  and  they  represent  no  obstacle  to 
ground-water  movement  between  Whirlwind  Valley  and  Sevier 
Desert.  However,  the  high  at  the  end  of  Profile  4  and  at  19  km 
on  Profile  5  (Figure  8)  may  be  evidence  of  a  bedrock  high  just 
east  of  Long  Ridge  and  Red  Knolls.  This  suggested  feature  could 
restrict  ground-water  migration  except  for  very  shallow  aqui¬ 


fers. 
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A1.0  GENERAL  PRINCIPLES  OF  THE  GRAVITY 
EXPLORATION  METHOD 

Al. 1  GENERAL 

A  gravity  survey  involves  measurement  of  differences  in  the 
gravitational  field  between  various  points  on  the  earth's 
surface.  The  gravitational  field  values  being  measured  are  the 
same  as  those  influencing  all  objects  on  the  surface  of  the 
earth.  They  are  generally  associated  with  the  force  which 
causes  a  1  gm  mass  to  be  accelerated  at  980  cm/sec^.  This 
force  is  normally  referred  to  as  a  1  g  force. 

Even  though  in  many  applications  the  gravitational  field  at  the 
earth's  surface  is  assumed  to  be  constant,  small  but  distin¬ 
guishable  differences  in  gravity  occur  from  point  to  point.  In 
a  gravity  survey,  the  variations  are  measured  in  terms  of 
milligals.  A  milligal  is  equal  to  0.001  cm/second2  or 
0.00000102  g.  The  differences  in  gravity  are  caused  by  geo¬ 
metrical  effects,  such  as  differences  in  elevation  and  latitude, 
and  by  lateral  variations  in  density  within  the  earth.  The 
lateral  density  variations  are  a  result  of  changes  in  geologic 
conditions.  For  measurements  at  the  surface  of  the  earth,  the 
largest  factor  influencing  the  pull  of  gravity  is  the  density  of 
all  materials  between  the  center  of  the  earth  and  the  point  of 
measurement . 


To  detect  changes  produced  by  differing  geological  conditions, 
it  is  necessary  to  detect  differences  in  the  gravitational  field 
as  small  as  a  few  milligals.  To  recognize  changes  due  to 
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geological  conditions,  the  measurements  are  "corrected"  to  ac¬ 
count  for  changes  due  to  differences  in  elevation  and  latitude. 

Given  this  background,  the  basic  concept  of  the  gravitational 
exploration  method,  the  anomaly,  can  be  introduced.  If,  instead 
of  being  an  oblate  spheroid  characterized  by  complex  density 
variations,  the  earth  were  made  up  of  concentric,  homogeneous 
shells,  the  gravi tational  field  would  be  the  same  at  all  points 
on  the  surface  of  the  earth.  The  complexities  in  the  earth's 
shape  and  material  distribution  are  the  reason  that  the  pull  of 
gravity  is  not  the  same  from  place  to  place.  A  difference  in 
gravity  between  two  points  which  is  not  caused  by  the  effects  of 
known  geometrical  differences,  such  as  in  elevation,  latitude, 
and  surrounding  terrain,  is  referred  to  as  an  "anomaly." 

An  anomaly  reflects  lateral  differences  in  material  densities. 
The  gravitational  attraction  is  smaller  at  a  place  underlain  by 
relatively  low  density  material  than  it  is  at  a  place  underlain 
by  a  relatively  high  density  material.  The  term  "negative 
gravity  anomaly"  describes  a  situation  in  which  the  pull  of 
gravity  within  a  prescribed  area  is  small  compared  to  the  area 
surrounding  it.  Low-density  alluvial  deposits  in  basins  such  as 
those  in  the  Nevada-Utah  region  produce  negative  gravity  anoma¬ 
lies  in  relation  to  the  gravity  values  in  the  surrounding 
mountains  which  are  formed  by  more  dense  rocks. 


The 

objective 

of  gravity 

exploration  is 

to  deduce  the 

variations 

i  n 

geologic 

conditions 

that  produce 

the  g  ravi ty 

a  noma  lies 

identified  during  a  gravity  survey. 
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Al. 2  INSTRUMENTS 

The  sensing  element  of  a  LaCoste  and  Romberg  gravimeter  is  a 
mass  suspended  by  a  zero-length  spring.  Deflections  of  the 
mass  from  a  null  position  are  proportional  to  changes  in  gravi¬ 
tational  attraction.  These  instruments  are  sealed  and  compen¬ 
sated  for  atmospheric  pressure  changes.  They  are  maintained  at 
a  constant  temperature  by  an  internal  heater  element  and  thermo¬ 
stat.  The  absolute  value  of  gravity  is  not  measured  directly  by 
a  gravimeter.  It  measures  relative  values  of  gravity  between 
one  point  and  the  next.  Gravitational  differences  as  small  as 
0.01  milligal  can  be  measured. 

Al. 3  FIELD  PROCEDURES 

The  gravimeter  readings  were  calibrated  in  terms  of  absolute 
gravity  by  taking  readings  twice  daily  at  nearby  USGS  gravity 
base  stations.  Gravimeter  readings  fluctuate  because  of  small 
time-related  deviations  due  to  the  effect  of  earth  tides  and 
instrument  drift.  Field  readings  were  corrected  to  account  for 
these  deviations.  The  magnitude  of  the  tidal  correction  was 
calculated  using  an  equation  suggested  by  Goguel  (1954): 

C  =  P  +  Ncos  0  (cos  0  +  sin  H  )  +  Scos  (cos  ff  -  sin  0) 
where  C  is  the  tidal  correction  factor,  P,  N,  and  S  are  time- 
related  variables,  and  #  is  the  latitude  of  the  observation 
point.  Tables  giving  the  'alues  of  P,  N,  and  S  are  published 
annually  by  the  European  Association  of  Exploration  Geophysi¬ 
cists. 
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The  meter  drift  correction  was  based  on  readings  taken  at  a 
designated  base  station  at  the  start  and  end  of  each  day.  Any 
difference  between  these  two  readings  after  they  were  corrected 
for  tidal  effects  was  considered  to  have  been  the  result  of 
instrumental  drift.  It  was  assumed  that  this  drift  occurred  at 
a  uniform  rate  between  the  two  readings.  Corrections  for  drift 
were  typically  only  a  few  hundredths  of  a  milligal.  Readings 
corrected  for  tidal  effects  and  instrumental  drift  represented 
the  observed  gravity  at  each  station.  The  observed  gravity 
values  represent  the  total  gravitational  pull  of  the  entire 
earth  at  the  measurement  stations. 

Al. 4  DATA  REDUCTION 

Several  corrections  or  reductions  are  made  to  the  observed 
gravity  to  isolate  the  portion  of  the  gravitational  pull  which 
is  due  to  the  crustal  and  near-surface  materials.  The  gravity 
remaining  after  these  reductions  is  called  the  “Bouguer 
Anomaly."  Bouguer  Anomaly  values  are  the  basis  for  geologic 
interpretation.  To  obtain  the  Bouguer  Anomaly,  the  observed 
gravity  is  adjusted  to  the  value  it  would  have  had  if  it  had 
been  measured  at  the  geoid,  a  theoretically  defined  surface 
which  approximates  the  surface  of  mean  sea  level.  The  dif¬ 
ference  between  the  "adjusted"  observed  gravity  and  the  gravity 
at  the  geoid  calculated  for  a  theoretically  homogeneous  earth  is 
the  Bouguer  Anomaly. 
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Four  separate  reductions,  to  account  for  four  geometrical 
effects,  are  made  to  the  observed  gravity  at  each  station  to 
arrive  at  its  Bouguer  Anomaly  value. 

a.  Free-Air  Effect:  Gravitational  attraction  varies  inversely 
as  the  square  of  the  distance  from  the  center  of  the  earth. 
Thus  corrections  must  be  applied  for  elevation.  Observed 
gravity  levels  are  corrected  for  elevation  using  the  normal 
vertical  gradient  of: 

FA  =  -0.09406  mg/ft  (-0.3086  mi  11 igals/meter ) 
where  FA  is  the  free-air  effect  (the  rate  of  change  of  gravity 
with  distance  from  the  center  of  the  earth).  The  free-air 
correction  is  positive  in  sign  since  the  correction  is  opposite 
the  effect. 

b.  Bouguer  Effect:  Like  the  free-air  effect,  the  Bouguer 

effect  is  a  function  of  the  elevation  of  the  station,  but  it 
considers  the  influence  of  a  slab  of  earth  materials  between 
the  observation  point  on  the  surface  of  the  earth  and  the 
corresponding  point  on  the  geoid  (sea  level) .  Normal  practice, 
which  is  to  assume  that  the  density  of  the  slab  is  2.67  grams 
per  cubic  centimeter  was  followed  in  these  studies.  The  Bouguer 
correction  (Bc),  which  is  opposite  in  sign  to  the  free-air 
correction,  was  defined  according  to  the  following  formula. 

Bc  =  0.01276  (2.67)  hf  (milligals  per  foot) 

Bc  =  0.04185  (2.61)  hm  (milligals  per  meter) 
where  hf  is  the  height  above  sea  level  in  feet  and  hm  is  the 
height  in  meters. 


' 
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c.  Latitude  Effect:  Points  at  different  latitudes  will  have 
different  "gravities"  for  two  reasons.  The  earth  (and  the 
geoid)  is  spheroidal,  or  flattened  at  the  poles.  Since  points 
at  higher  latitudes  are  closer  to  the  center  of  the  earth  than 
points  near  the  equator,  the  gravity  at  the  higher  latitudes  is 
larger.  As  the  earth  spins,  the  centrifugal  acceleration 
causes  a  slight  decrease  in  gravity.  At  the  higher  latitudes 
where  the  earth's  radii  are  smaller,  the  centrifugal  accel¬ 
eration  diminishes.  The  gravity  formula  for  the  Geodetic 
Reference  System,  1967,  gives  the  theoretical  value  of  gravity 
at  the  geoid  as  a  function  of  latitude.  It  is: 

g  =  978.0381  (1  +  0.0053204  sin2  jrf  -  0.0000058  sin22(rf)  gals 
where  g  is  the  theoretical  acceleration  of  gravity  and  is 
the  latitude  in  degrees.  The  positive  term  accounts  for  the 
spheroidal  shape  of  the  earth.  The  negative  term  adjusts  for 
the  centrifugal  acceleration. 

The  previous  two  corrections  (free  air  and  Bouguer)  have  ad¬ 
justed  the  observed  gravity  to  the  value  it  would  have  had  at 
the  geoid  (sea  level) .  The  theoretical  value  at  the  geoid  for 
the  latitude  of  the  station  is  then  subtracted  from  the  adjusted 
observed  gravity.  The  remainder  is  called  the  Simple  Bouguer 
Anomaly  (SBA) .  Most  of  this  gravity  represents  the  effect  of 
material  beneath  the  station,  but  part  of  it  may  be  due  to 
irregularities  in  terrain  (upper  part  of  the  Bouguer  slab)  away 


from  the  station. 


d.  Terrain  Effect;  Topographic  relief  around  the  station  has 
a  negative  effect  on  the  gravi tational  force  at  the  station.  A 
nearby  hill  has  upward  gravitational  pull  and  a  nearby  valley 
contributes  less  downward  attraction  than  a  nearby  material 
would  have.  Therefore,  the  corrections  are  always  positive. 
Corrections  are  made  to  the  SBA  when  the  terrain  effects  were 
0.1  milligal  or  larger.  Terrain  corrected  Bouguer  values  are 
called  the  Complete  Bouguer  Anomaly  (CBA).  When  the  C BA  is 
obtained,  the  reduction  of  gravity  at  individual  measurement 
points  (stations)  is  complete. 

A1.5  INTERPRETATION 

The  first  step  in  interpretation  is  to  separate  the  portion  of 
the  CBA  that  might  be  caused  by  the  lightweight,  basin-fill 
material  overlying  the  heavier  bedrock  material  which  forms  the 
surrounding  mountains  and  presumably  the  basin  floor.  Since  the 
valley-fill  sediments  are  absent  at  the  stations  read  in  the 
mountains,  the  CBA  values  at  these  bedrock  stations  are  used  as 
the  basis  for  constructing  a  regional  field  over  the  valley.  A 
regional  field  is  an  estimation  of  the  values  the  CBA  would  have 
had  if  the  light  weight  sediments  (the  anomaly)  had  not  been 
there . 

The  difference  between  the  CBA  and  the  regional  field  is 
called  the  "residual"  field  or  residual  anomaly.  The  residual 
field  is  the  interpreter's  estimation  of  the  gravitational 
effect  of  the  geologic  anomaly.  The  zero  value  of  the  residual 
anomaly  is  not  exactly  at  the  rock  outcrop  line  but  at  some 
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distance  on  the  "rock"  side  of  the  contact.  The  reason  for  this 
is  found  in  the  explanation  of  the  terrain  effect.  There  is  a 
component  of  gravitational  attraction  from  material  which  is 
not  directly  beneath  a  point. 

If  the  "regional"  is  well  chosen,  the  magnitude  of  the  residual 
anomaly  is  a  function  of  the  thickness  of  the  anomalous  (fill) 
material  and  the  density  contrast.  The  density  contrast  is  the 
difference  in  density  between  the  alluvial  and  bedrock  material. 
If  this  contrast  were  known,  an  accurate  calculation  of  the 
thickness  could  be  made.  In  most  cases,  the  densities  are  not 
well  known  and  they  also  vary  within  the  study  area.  In  these 
cases,  it  is  necessary  to  use  typical  densities  for  materials 
similar  to  those  in  the  study  area. 

If  the  selected  average  density  contrast  is  smaller  than  the 
actual  density  contrast,  the  computed  depth  to  bedrock  will  be 
greater  than  the  actual  depth  and  vice-versa.  The  computed 
depth  is  inversely  proportional  to  the  density  contrast.  A  ten 
percent  error  in  density  contrast  produces  a  ten  percent  error 
in  computed  depth.  An  iterative  computer  program  is  used  to 
calculate  a  subsurface  model  which  will  yield  a  gravitational 
field  to  match  (approximately)  the  residual  gravity  anomaly. 
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